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PREFACE 


The  Dynamic  Environment  Simulator  (DES)  is  a  three  axes  centrifuge  device  which  has  been 
used  to  replicate  agile  flight  or  super  maneuvers  that  an  aircraft  may  fly  in  a  combat  scenario. 
In  order  to  simulate  these  maneuvers,  it  requires  the  centrifuge  to  utilize  simultaneous 
multiaxes  motion.  In  this  device,  a  pilot  will  sit  on  the  seat  in  the  third  and  innermost 
gimbal  (the  cab)  and  be  subjected  to  simultaneous  three  axes  motion.  As  all  three  axes 
rotate,  the  subject  will  not  only  experience  the  motion  programmed  to  emulate  the 
supermaneuver,  but  also  experience  motion  fields  generated  by  Coriolis  components.  These 
Coriolis  components  produce  unusual  and  unnatural  accelerations  (which  induce  forces)  and 
reduce  the  credibility  of  the  motion  simulation  to  emulate  the  true  flight  scenario. 

Coriolis  accelerations  are  a  natural  property  arising  from  physics  and  dynamics  and  occur 
whenever  a  coordinate  frame  moves  relative  to  another  moving  coordinate  frame.  The  DES 
is  a  physical  system  and,  as  such,  each  gimbal  is  connected  to  another  supporting  gimbal. 
Thus  the  very  nature  of  the  physical  construction  of  the  DES  provides  a  framework  to 
produce  Coriolis  type  accelerations. 

In  order  to  work  on  a  means  of  reducing  the  Coriolis  action  produced  by  the  centrifuge 
motion  simulator,  the  first  step  is  to  accurately  define  the  equations  of  motion  of  the  DES. 
Once  these  equations  of  motion  are  obtained,  one  can  then  extract  those  terms  that  contribute 
only  to  the  Coriolis  component.  Knowing  these  terms,  the  next  step  is  to  minimize  their 
magnitude  by  inserting  into  the  commands  of  the  actuators  an  additional  signal  to  drive  the 
centrifuge  simulator  in  such  a  manner  as  to  mitigate  these  untoward  effects.  This 
minimization  of  Coriolis  components  also  has  a  price  in  that  it  compromises  the  other 
commands  inserted  into  the  actuators,  thereby  affecting  the  total  motion  field  produced. 

Thus  it  is  necessary  to  both  characterize  the  extent  of  the  motion  field  action  caused  by  the 
Coriolis  accelerations  as  well  as  when  this  action  is  minimized.  Consequently,  it  is  also 
necessary  to  quantify  the  impact  of  the  minimization  of  the  Coriolis  on  the  overall  credibility 
of  the  motion  simulation. 


To  minimize  the  Coriolis  effects,  an  optimization  procedure  is  given  in  this  report  which 
allows  the  Control  Engineer  the  opportunity  to  insert  appropriate  input  commands  into  the 
actuators  of  the  centrifuge  simulator.  These  input  commands  turn  out  to  have  a  special 
feedback  nature  which  not  only  provides  both  a  method  to  accurately  control  a  mechanical 
system  such  as  the  DES,  but  also  allows  the  monitoring  and  regulation  of  the  response  of  this 
motion  simulator. 
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INTRODUCTION 


The  DBS  (Dynamic  Environment  Simulator)  is  a  three  axes  centrifuge  simulator  which  can 
be  used  to  emulate  motion  fields  as  experienced  in  agile  or  supermaneuverable  aircraft  [1,2]. 
Figure  1  illustrates  a  schematic  diagram  of  this  mechanical  system.  In  order  to  generate 
complex  motions  that  appear  in  such  flight  scenarios,  one  or  more  axes  of  the  DBS  must  be 
operated  simultaneously  to  develop  the  desired  motion  field  profile  to  appear  at  the  seat 
located  within  the  cab  depicted  in  Figure  1. 
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Figure  1 .  A  Schematic  Diagram  of  the  DBS 

When  a  coordinate  frame  moves  relative  to  another  moving  frame  (which  occurs  in  the  DBS 
since  the  physical  gimbals  are  interconnected  to  one  another),  this  may  introduce  Coriolis 
acceleration  components.  What  is  unusual  about  the  Coriolis  induced  components  are  that 
they  act  in  a  direction  perpendicular  to  normal  induced  accelerations  and  motion.  When 
humans  are  exposed  to  such  phenomena,  they  feel  these  motions  are  very  unrealistic  since 
the  forces  they  experience  are  orthogonal  to  the  direction  of  movement.  Since  the  goal  in 
motion  simulation  is  to  generate  both  linear  and  angular  accelerations  that  feel  realistic,  the 
Coriolis  component  is  detrimental  to  the  production  of  good  motion  in  a  perceptual  sense. 

The  purpose  of  this  study  is  to  first  quantify  the  Coriolis  component  of  the  motion  produced 
in  the  DBS.  Secondly,  it  is  desired  to  work  on  a  means  of  reducing  this  effect  without  too 
much  compromise  of  the  other  motions  that  are  chosen  to  be  emulated.  First  an  argument  is 
made  on  the  existence  of  such  motions  as  Coriolis  and  some  examples  are  given  to  illustrate 
this  natural  artifact  of  relative  motion.  These  motions  are  then  quantified,  and  then  a  means 
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of  reducing  such  motions  is  suggested.  A  computer  simulation  demonstrates  the  extent  of  the 
reduction  of  the  Coriolis  component  by  the  algorithms  employed  herein.  We  first  discuss  the 
existence  of  such  untoward  motions  and  how  they  are  a  natural  product  of  the  physics  and 
dynamics  of  mechanical  systems  such  as  the  DES. 


Definition  of  Coriolis 


The  classical  definition  of  Coriolis  acceleration  is  well-known  since  its  early  definition  by  the 
French  mathematician,  G.  Coriolis.  The  first  example  in  this  report  is  used  to  define 
Coriolis  in  a  classical  sense. 


Example  1  -  The  Classical  Definition  of  Coriolis 

If  a  bar  rotates  (Figure  2)  with  an  angular  velocity  W  which  has  a  vector  direction  out  of  the 
paper  and  a  slider  moves  along  the  arm  with  linear  velocity  V,  the  action  of  the  Coriolis 
acceleration  a(corioiis)  which  is  induced  on  the  slider  satisfies  the  vector  cross  product 
relationship: 


a(CorioUs)  =  2  W  X  V 


(1) 


where  a(corioiis)  is  in  a  direction  orthogonal  to  both  the  arm  and  the  W  vector  as  indicated  in 
Figure  2.  If  a  human  subject  were  positioned  on  the  slider  on  the  rotating  arm,  and  were  to 
make  any  body  motion  in  the  direction  of  the  slider  movement  (V),  an  acceleration  a^conoiis) 
would  act  on  the  subject  in  the  direction  of  a(Corioiis)  as  depicted  in  Figure  2.  This  is  very 
counterintuitive  for  a  human  subject  to  experience  an  acceleration  (or  resulting  force)  on  his 
body  acting  in  this  direction,  orthogonal  to  the  moving  slider. 
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Figure  2.  The  Classical  Definition  of  Coriolis 
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More  specifically,  if  a  human  subject  were  positioned  on  the  arm  with  his  head  closest  to  the 
center  of  rotation,  and  if  he  were  to  extend  his  arm  radially  outward,  a  force  would  be  felt  to 
his  left  side.  If  he  reverses  this  motion,  he  feels  a  force  to  the  right.  This  is  also  very 
counterintuitive  because  the  directions  change.  For  further  discussion,  use  will  be  made  of 
Figure  3  which  illustrates  the  body  centered  coordinate  system  notation  normally  applied  in 
sustained  acceleration  studies  and  will  be  adopted  as  the  reference  frame  in  this  report. 
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Figure  3.  The  Body  Axis  Coordinate  System 


Figure  4.  The  Two  Disk  Example  to  Mathematically  Derive  Coriolis 
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A  second  example  is  now  given  which  illustrates  a  mathematical  basis  why  Coriolis 
accelerations  always  appear  on  systems  where  a  frame  experiences  movement  relative  to 
another  ambulatory  frame  (usually  inertial). 

Example  2  -  The  Mathematical  Basis  of  Coriolis  Motion 

Figure  4  depicts  an  example  from  dynamics  which  provides  a  clear  mathematical 
understanding  of  the  origin  of  Coriolis  accelerations  and  their  subsequent  forces.  In  Figure 
4,  the  large  disk  (Disk  1)  of  radius  R1  rotates  at  angular  velocity  Wl.  The  smaller  disk 
(Disk  2)  of  radius  r2  has  the  same  center  and  rotates  at  angular  velocity  W2  relative  to  Disk 
1 .  It  is  desired  to  calculate  the  accelerations  appearing  at  the  point  A  on  the  smaller  disk. 
The  centripetal  acceleration  acting  at  a  point  moving  in  a  circle  is  given  by: 

acceleration(,e„,ripetoi)  =  Radius  *  (2) 

For  simplicity,  this  example  will  assume  Wl  and  W2  are  constant,  thus  eliminating  any 
tangential  component  of  acceleration  appearing.  Hence,  the  centripetal  acceleration  is  the 
main  acceleration  of  interest  at  the  point  A.  Since  both  disks  rotate  about  the  same  center 
and  the  smaller  disk  rotates  on  top  of  the  larger  disk,  the  total  angular  velocity  of  point  A 
with  respect  to  the  inertial  frame  is  Wl  +  W2.  There  are  two  ways  to  attempt  to  calculate 
the  net  radial  acceleration  appearing  at  point  A. 

Method  1  -  Attempt  to  Add  the  Accelerations  Linearly 

At  first  glance,  one  may  try  to  add  the  two  accelerations  separately  in  a  linear  manner  by 
taking  into  account  the  movement  of  both  bodies.  Proceeding  in  this  manner,  the 
acceleration  induced  on  point  A  by  the  larger  disk  would  be  given  by: 

centripetal  acceleration  due  to  Disk  1  =  r2  Wl^  (3) 

since  the  point  A  is  r2  units  from  the  center  of  rotation  of  disk  1.  Next,  one  would  calculate 
the  acceleration  induced  at  point  A  due  to  the  rotation  of  the  smaller  disk. 

centripetal  acceleration  due  to  Disk  2  =  r2  W2^  (4) 

Thus,  if  this  procedure  is  followed,  the  total  acceleration  acting  on  point  A  would  be 
erroneously  calculated  as: 

total  acceleration  on  point  A  =  r2  Wl^  +  r2  W2^  (5) 

This  is  not  correct  because  the  Coriolis  effects  have  been  neglected.  The  missing  Coriolis 
term  is: 

missing  Coriolis  term  =  2  *  Wl  r2  W2  (6) 
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and  acts  in  a  radial  direction  on  the  smaller  disk  at  point  A.  To  see  how  this  missing  term 
arises,  the  correct  solution  to  this  problem  is  now  presented  to  illustrate  its  origin. 

Method  2  -  Calculation  of  Accelerations  from  the  Inertial  Frame 

If  one  now  views  the  system  in  Figure  4  from  an  inertial  frame  perspective,  the  point  A  has 
a  total  angular  velocity  of  W1  +  W2.  Thus,  the  calculation  of  the  total  radial  acceleration 
acting  on  the  point  A  is  given  by: 

total  centripetal  acceleration  of  point  A  =  (Wl+W2)^*r2  (7) 

This  differs  from  the  solution  given  in  equation  (5)  but  is  the  correct  result.  With  some 
study,  the  missing  term  can  be  attributed  to  the  fact  that: 

WP  +  W22  is  not  equal  to  (W1  +  'Wlf  (8) 

Thus,  Coriolis  accounts  for  the  cross  term  which  is  missing  when  the  square  of  the  sum 
differs  from  the  sum  of  the  squares.  This  difference,  however,  may  comprise  a  very 
significant  part  of  the  calculation  of  the  total  acceleration  that  appears  at  point  A. 

There  are  some  situations  in  which  it  is  physically  possible  to  experience  the  untoward 
artifact  of  Coriolis  accelerations.  Figure  5  illustrates  a  simple  playground  device  for  children 
which  enables  anyone  to  perceive  this  variable  and  validate  its  existence.  If  the  device 
rotates  counter  clockwise  and  a  subject  sits  (looking  inward)  on  the  rim  of  the  apparatus 
making  an  arm  or  foot  movement  radially  inward  will  result  in  an  experience  of  an 
acceleration  (and  force)  to  the  subject’s  left.  Conversely,  making  an  arm  or  foot  movement 
radially  outward  will  result  in  an  experience  of  an  acceleration  (and  force)  to  the  right.  If 
the  direction  of  rotation  of  the  playground  device  is  reversed,  opposite  accelerations  and 
forces  are  perceived. 

Finally,  the  use  of  the  term  "Coriolis"  should  be  distinguished  from  its  common  usage 
occurring  in  studies  on  spatial  orientation.  This  has  a  different  meaning  and  refers  to  an 
illusion  or  misperception  of  attitude  information.  For  example,  from  a  subject’s  perspective, 
if  he  moves  his  head  out  of  one  plane  of  rotation  and  into  another,  as  can  occur  on  a  human 
centrifuge,  the  semi-circular  canals  of  the  human  vestibular  system  can  be  excited  in  such  a 
manner  as  to  elicit  conflicting  motion  cues.  The  "Coriolis  illusion,"  as  it  is  called,  is 
reduced  in  pilots  by  having  them  restrict  head  motions  while  their  aircraft  is  turning  or 
accelerating.  Reference  [3]  addresses  the  detection  of  such  illusory  phenomenon. 
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Figure  5.  A  Physical  System  To  Experience  Coriolis  Acceleration 

The  Coriolis  artifact  just  described  is  significantly  larger  in  a  motion  simulator  as  compared 
to  the  actual  aircraft  because  the  turn  radius  difference  is  two  orders  of  magnitude  larger  in 
the  aircraft  as  compared  to  the  simulator.  This  type  of  motion  mismatch  degrades  the 
fidelity  of  the  motion  simulator  considered  herein  and  limits  its  use  to  very  restricted  motion 
studies.  It  is  desired  in  this  report  to  modify  the  joint  commands  into  the  actuators  of  the 
simulator  to  eliminate,  or  at  least  reduce,  some  of  these  Coriolis  accelerations  which  are 
produced  in  the  ground  base  motion  simulator.  In  this  manner  the  motion  simulator  will 
produce  a  motion  field  which  is  more  credible  to  human  subjects,  yet  utilize  the  full 
capability  of  the  multi  axes  system  of  the  centrifuge  motion  simulator. 

The  next  step  is  to  describe  the  DES  in  more  explicit  detail  and  then  derive  the  equations  of 
motion. 


The  Dynamic  Environment  Simulator  (DES) 

The  approach  taken  here  will  be  to  model  this  mechanical  system  within  the  context  of  a 
robotic  system.  This  is  due  to  many  new  developments  in  the  study  of  robotic  manipulators 
that  will  help  in  the  determination  of  the  Coriolis  terms  that  appear  in  the  equations  of 
motion. 
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If  one  wishes  to  model  a  three  axes  centrifuge  simulator  as  a  robotic  system  (Figure  1),  one 
approach  is  to  assume  the  centrifuge  is  a  multiple  link  system  with  three  revolute  axes  in 
which  each  axis  rotates  360  degrees.  The  human  (a  pilot  involved  in  an  experiment  of 
unusual  acceleration  stress  or  rotations)  sits  at  the  end  effector  (seat  in  the  cab)  and  his  body 
is  subjected  to  the  induced  motion  field.  When  modeling  a  centrifuge  simulator  as  a  robotic 
system,  the  main  distinction  between  this  problem  and  a  traditional  robotics  problem  involves 
the  size  of  the  mechanical  system  which,  in  this  case,  is  prohibitively  large  (over  24  feet  in 
height  and  over  180  tons  in  weight)  with  a  work  volume  exceeding  9x10^  feet^. 

Once  this  problem  is  defined  within  the  context  of  a  robotic  system,  the  forward  and  inverse 
kinematics  problems  which  help  define  the  equations  of  motion  can  be  obtained  using 
standard  Denavit-Hartenberg  formulation.  As  a  first  step,  it  is  required  to  classify  all 
centrifuge  parameters  (both  fixed  and  rotating  quantities)  that  will  influence  the  equations  of 
motion,  and  thus  the  calculated  Coriolis  induced  accelerations. 

Table  1  illustrates  the  physical  parameters  of  the  centrifuge  simulator  depicted  in  Figure  1. 
The  last  six  fixed  variables  (L„Ly,L„L„  ©,y,  and  refer  to  the  seat  in  the  cab  displayed 
in  Figures  1  and  3. 


Table  1.  Fixed  and  Rotational  Parameters  of  the  DES 


Fixed  Rotational 


Variable 

Dimensions 

Variable 

Max 

Min 

Dimensions 

L. 

12.27  Feet 

9' 

2  Pi 

-2  Pi 

rad. 

4 

19.0  Feet 

5.86 

-5.86 

rad/s 

4 

3.90  Feet 

02 

0.63 

-0.63 

rad/s/s 

Ly 

0.012  Feet 

©2 

2  Pi 

-2.0  Pi 

rad 

L, 

0.012  Feet 

3.14 

-  Pi 

rad/s 

L. 

2.76  Feet 

e. 

2.5 

-2.5 

rad/s/s 

©.y 

0<  ©_y<2  Pi  rad 

©3 

2  Pi 

-2  Pi 

rad 

0<e„<Pi/2  rad 

93 

3.14 

-3.14 

rads/s 

©1 

5.0 

-5.0 

rad/s/s 

where  Pi=3.1416,  ©i  (arm),  ©2  (fork),  and  ©3  (cab)  are  the  three  revolute  joints.  The 
variables  Q^y  and  refer  to  rotations  of  the  pilot’s  seat  (yaw  and  pitch)  at  the  end  effector. 
The  variables  L,  and  L2  refer  to  centrifuge  link  lengths;  the  variables  L,,  Ly,  L^,  and  Lt  refer 
to  link  lengths  associated  with  the  pilot’s  seat  at  the  end  effector  (cf.  Figures  1,3). 

From  a  pilot’s  perspective,  if  the  seat  is  configured  in  the  tangential  direction  (©sy  =  Pi/2 
radians),  the  traditional  roll,  pitch,  and  yaw  angles  would  be  specified  by:  roll  =  ©3,  pitch 
=  -  ©2,  and  yaw  =  ©^y.  It  is  also  important,  in  describing  the  DES  as  a  robotic 
manipulator,  to  clarify  how  the  controlling  action  of  the  different  axes  is  achieved. 
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The  Control  Characteristics  of  the  Device 


The  first  axis  in  Figure  1  (the  arm  described  by  the  axis  ©i  or  Wl)  is  a  velocity  control 
system  driven  by  three  electric  motors  physically  connected  in  tandem.  Thus,  velocity 
command  inputs  are  compared  to  the  angular  velocity  (tachometer)  measurement  of  the  actual 
arm  speed  producing  a  velocity  error.  Acting  in  a  feedback  manner,  this  velocity  error  is 
multiplied  by  a  gain  and  this  signal  augments  the  torque  produced  for  that  axis  to  drive  the 
attained  velocity  to  the  command  signal.  The  second  axis  (the  fork  with  axis  ©2)  is  actuated 
by  a  separate  DC  electric  motor  but  is  a  position  control  system.  In  this  case,  the  input 
consists  of  position  command  signals  (voltages)  and  position  measurements  are  obtained  from 
a  pair  of  resolvers,  thus  delivering  a  position  error  for  use  in  a  controller  design.  The  third 
axis  (the  cab  with  axis  ©3)  is  also  a  position  control  system,  but  is  actuated  by  a  hydraulic 
motor  and  has  similar  operation  to  the  fork  control  system. 

In  order  to  specify  the  actual  motion  produced  by  the  flight  simulator,  the  kinematics  aspects 
of  this  manipulator  are  now  described  in  more  specific  detail. 


The  Kinematic  Description  of  the  DBS 

In  [1]  a  description  of  the  forward  kinematics  of  this  mechanical  system  was  derived  using 
the  frame  assignments  (Denevit-Hartenberg)  from  Craig  [4].  This  detailed  derivation 
involving  a  total  of  17  coordinate  systems  was  used  to  display  many  different  possible  seat 
orientations  and  translations  as  illustrated  in  Table  2.  The  purpose  of  using  such  a  large 
number  of  coordinate  frames  in  this  analysis  of  the  manipulator  system  was  to  have 
additional  freedom  to  investigate  possible  optimization  of  the  system’s  configuration 
parameters. 

Other  kinematic  derivations  are  much  simpler,  e.g.  [5]  and  are  useful  for  other  situations. 
We  display  here  the  parameters  obtained  by  viewing  the  DES  within  the  context  of  a  robotic 
manipulator  system  using  the  full  Denevit-Hartenberg  formulation. 
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Table  2.  The  Denavit-Hartenberg  Parameters  For  The  DES 


From 

Distance  Along 

Distance  Along 

Angle  Fr 

Angle 

zi-1  to  zi 

zi-1  to  zi 

Xi-1  to  Xi 

Xi-1  to  > 

Link 

about  Xi-1 

along  Xi-1 

along  Zi 

about  Zi 

i 

=  ai-1 

=  Ai-1 

=  di 

=  ©i 

1 

0 

0 

0 

-©1 

2 

0 

L2 

0 

0 

3 

0 

0 

0 

Pi/2 

4 

Pi/2 

0 

0 

0 

5 

0 

0 

0 

©2 

6 

0 

0 

0 

Pi/2 

7 

-Pi/2 

0 

0 

0 

8 

0 

0 

0 

©3 

9 

0 

-Lx 

+Lz 

0 

10 

0 

0 

0 

-  Pi/2 

11 

0 

Ly 

0 

0 

12 

Pi/2 

0 

0 

0 

13 

0 

0 

0 

“  ^sy 

14 

-Pi/2 

0 

0 

0 

15 

0 

0 

0 

-  ©St 

16 

Pi/2 

0 

0 

0 

17 

0 

0 

-Lt 

0 

The  transformation  matrix  Tq*’  relating  the  inertial  frame  to  the  seventeenth  coordinate 
system  can  be  derived  (cf.  [6]): 


T„, 

Tn, 

Tb, 

Th 

T21, 

'^22i 

^23* 

T24 

T3,, 

T32> 

T33> 

T34 

(9) 

T41, 

T42> 

T43> 

T44 

where  Tq*’  is  the  Denevit-Hartenberg  transformation  used,  e.g.  in  [4],  with  relevant 
components: 


Til  —  si  c2  S(sy)  C(s,)  +  si  s2  c3  S(s() 

-si  s2  s3  C(,y)  C(,,)  -1-  cl  s3  S(3i)  -1-  clc3C(,y)C(„) 

(10a) 

Ti2  =  -si  c2  C(sy)  -si  s2  s3  S(,y)  -f  cl  c3  S(,y) 

Ti3  -sl  c2  s^jy)  S(5t)  -f"  si  s2  c3  c^^jj 

(10b) 

-l-  sl  s2  s3  C(5y)  S(5()  -h  cl  s3  C(g{)-clc3C(sy)S(50 

(10c) 

Ti4  =  sl  c2  Lt  S(sy)  S(st)  -Lz  sl  c2-sl  s2  c3  Lt  C(si) 

-sl  s2  s3  Lt  C(sy)  S(st)  -sl  s2  Ly  s3  -cls3LtC(s,) 

+  cl  Lx  s3  +cl  Lt  C(sy)  S(st)  c3+  cl  Lyc3+cl  L2+Lxc3sls2  (lOd) 
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T21  —  cl  c2  S(sy)  C(gt)  "l“cl  s2  c3 

“Cl  S2  S3  ^(st) 

(lOe) 

T22  =  -  cl  c2  C(sy)  -  cl  s2  s3  S(gy)  -si  c3  s^gy) 

^23  ~  -  cl  c2  S(sy)  S^gt)  "1"  cl  s2  c3  C(gQ 

(lOf) 

-|-  cl  S2  S3  C^sy)  ^(st)  ^(st)  sl^3c^gy^S(g^^ 

(lOg) 

T24  =  cl  c2  Lt  S(gy)  S(gt)  -  cl  c2  Lz  -cls2c3  LtC(gt) 

+  cls2c3  Lx  -cls2  Lt  C(sy)S(gt)S3  -  cls2  Lys3 
+  si  s3  Lt  C(gt)  -  si  Lx  s3  -si  Lt  S(gt)C3C(gy) 

-Ly  si  c3  -L2  si 

(lOh) 

^31  “  s2  S(sy)  C(st)  -  c2  c3  S(gj)  +  c2  s3C(sy)C(st) 

(lOi) 

T32  ^  ’^2  C(sy)  +  C2  S3  S(sy) 

(lOj) 

T33  =  “S2  S(sy)  S(st)  -c2  c3  C(st)  -  c2  s3C(sy)S(st) 

(10k) 

X34  ~  s2  Lt  s^gy^  s^gj^  ”S2  Lz  “I"c2  c3  Lt  c^gj^ 

-•c2  c3  Lx  +  c2  s3  Lt  C(sy)  S(st)  +  c2  Ly  s3 

(lOL) 

where  s  and  c  refer  to  the  standard  sine  and  cosine  notation  and  for  the  subscript  notation. 

those  terms  with  s  refer  only  to  the  seat  variables. 

To  complete  the  kinematic  relationships,  the  relevant  details  used  in  the  derivation  can  be 
briefly  summarized  as  follows: 


Forward  Kinematic  Relationships 

Let  X  be  a  3x1  position  vector  in  Cartesian  space  and  ©  a  3x1  joint  space  vector.  The 
forward  kinematics  relationship  between  X  and  0  can  be  specified  via: 

X  =  f(  0  )  (11) 

where  the  function  f(.)  [6]  for  a  centrifuge  motion  simulator  is  an  extremely  complicated 
function  of  the  3  joint  variables  ©j,  (i= 1,3)  and  their  derivatives.  If  X  is  a  linear  position 
vector  and  denoting  its  3x1  dimensional  velocity  vector  as  V,  then: 

V  =  (d/dt)  X  =  J1  (d/dt)  ©  (12) 

as  the  linear  velocity  relationship  and  the  dot  or  derivative  indicates  time  differentiation.  J1 
is  the  3x3  manipulator  Jacobian  matrix  associated  with  linear  velocities. 
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If  the  vector  W  represents  the  angular  velocity  of  a  coordinate  frame  located  at  the  end 
effector  within  the  simulator,  then  the  representation: 

W  =  J2  (d/dt)  ©  (13) 

is  appropriate  to  describe  the  rotational  velocities  in  Cartesian  space.  Here  J2  is  another  3x3 
manipulator  Jacobian  matrix  associated  with  the  rotational  velocities.  Combining  equations 
(12,13),  it  is  seen  that: 

H  =  [n]  (d/dt)  0  =  J  (d/dt)  0  (14) 

where  J  is  defined  as  the  overall  6x3  manipulator  Jacobian  for  this  system.  J  can  be 
calculated  in  a  variety  of  ways  (Orin  and  Schrader,  [7]). 

This  report  addresses  the  application  of  the  motion  simulator  to  disorientation  studies  where 
the  principal  focus  involves  the  orientation  Jacobian  J2.  Figure  3  illustrates  the  body 
centered  coordinate  system  located  at  the  end  effector  in  which  all  the  final  motion 
calculations,  involving  rotations  and  linear  accelerations,  are  performed. 

In  the  work  considered  herein,  two  minimum  Coriolis  inverse  kinematics  algorithms  are 
considered  to  drive  the  motion  simulator.  The  first  algorithm  contains  zero  Coriolis;  the 
second  algorithm  minimizes  the  Coriolis.  The  purpose  of  the  development  of  these 
algorithms  is  to  modify  the  input  kinematic  commands  of  a  motion  simulator,  similar  to  a 
centrifuge,  to  make  the  resulting  motion  fields  more  realistic  to  humans. 


Calculation  of  Coriolis  Induced  on  the  DES 

As  mentioned  previously,  let  W(t)  represent  the  total  angular  velocity  of  the  third,  fixed  body 
(the  cab)  desired  at  the  end  point  of  the  motion  simulator  in  Figure  1.  This  report  addresses 
just  the  orientation  problem  to  illustrate  the  procedure.  For  a  three  axes  centrifuge  motion 
simulator  (modeled  as  a  robot),  one  could  express  the  3x1  vector  (d/dt)  W(t)  (the  total  vector 
of  angular  accelerations)  measured  in  the  same  coordinate  frame  as  W(t)  in  terms  of  the  unit 
vectors  of  the  inertial  frame  (ioJojko)  in  the  following  manner: 

(d/dt)W(t)  =  [A3(t)io,  B3(t)jo,  C3(t)kof  (15) 

where  it  can  be  shown  that  for  the  centrifuge  simulator: 

A3(t)  =  (dW)©2  cl  -  (d/dt) ©2  [(d/dt)©,]  si 
-(dW)©3sl  c2  -  (d/dt) ©3  [(d/dt)©,]  cl  c2 

+  (d/dt)©3  [(d/dt)©2]  si  s2  (16) 
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B3(t)  =  -  (dW)Q2  si  -  (d/dt)©2  [(d/dOGJ  cl 
-  (dW)e3  cl  c2  +  (d/dOGj  [(d/dOGj  si  c2 

+  (d/dOGj  [(d/dt)G2]  cl  s2  (17) 

C3(t)  =  -  (dW)G,  -  (dW)G3  s2 -(d/dt)G3  [(d/dt)G2]  c2  (18) 

where  i= 1,3  is  a  vector  describing  the  joint  angles  of  the  centrifuge  simulator  and  si 
and  ci  refer  to  sine  and  cosine,  respectively.  The  orientation  Jacobian  matrix  J2  was  defined 
in  equation  (13)  and  the  inverse  kinematics  relationship  for  orientation  commands  by  which 
the  joint  angles  can  be  integrated  forward  (since  W(t)  is  specified  a  priori)  is  given  by: 

(d/dt)  G  =  [J2]-'  W  (19) 

which  is  the  inverse  kinematics  solution  of  interest  in  this  study.  The  assumption  in  this 
report  is  that  the  orientation  Jacobian  J2  remains  in  a  reasonable  full  rank  configuration.  In 


[1]  this  problem  was  addressed  where  this  assumption  was  relaxed. 

It  is  noted  that  (d/dt)  W(t)  can  be  broken  up  into  two  terms: 

(d/dt)W(t)  =  (d/dt)W„,(t)  +  (d/dt)We(t)  (20) 

where  the  subscripts  nc  and  c  refer  to  noncoriolis  and  Coriolis  terms,  respectively.  For  the 
motion  simulation  problem  this  can  be  more  explicitly  written  in  the  inertial  coordinate 
frame: 

(d/dt)W,(t)  =  [AB3(t)io,  BB3(t)jo,  CB3(t)ko]T  (21) 

and  the  expressions  for  AB3,  BB3,  and  CB3  can  be  easily  obtained  from  equations  (16-18): 

AB3(t)  =  -(d\dt)G2  [(d\dt)G,]  si  -  (d\dt)G3  [(d\dt)G,]  clc2 

-l-(d\dt)G3  [(d\dt)G2]  si  s2  (22) 

BB3(t)  =  -(d\dt)G2  [(d\dt)G,]  cl  -I-  (d\dt)G3  [(d\dt)G,]  slc2 

+  (d\dt)  G 3  [(d\dt)  G 2]  cl  s2  (23) 

CB3(t)  =  -  (d\dt)G3  [(d\dt)G2]  c2  (24) 


The  efforts  in  this  report  will  focus  on  equation  (20)  dealing  with  orientation  variables 
because  Coriolis  is  so  easily  quantified.  This  term  is  what  predominately  affects  pilots  in  a 
motion  simulator.  It  is  now  possible  (assuming  a  finite  time  t  in  [to,  tj  )  to  define  two 
optimization  problems  related  to  minimizing  Coriolis  accelerations  (and  their  resulting  forces) 
appearing  at  the  endpoint  of  the  motion  simulator.  When  viewed  within  the  context  of  an 
optimization  control  problem  (Lee  and  Markus,  [8],  and  Bryson  and  Ho,  [9]),  these  problems 
are  linear-quadratic  in  nature.  In  the  next  section,  Gj^j  performs  the  role  of  state  vector 
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which  is  the  adapted  joint  commands  and  p(t)  is  the  control  variable  which  needs  to  be 
determined.  First  the  problem  will  be  formulated  in  attitude  space  with  a  hard  constraint  of 
zero  Coriolis  to  appear  at  the  end  effector. 


A  Means  to  Reduce  Coriolis  on  the  DES 
Method  1  -  Zero  Coriolis  in  an  Attitude  Space 

The  first  attempt  at  this  problem  is  to  find  a  time  history  ©^(t)  which,  for  a  robotic 
manipulator  system,  becomes  the  modified  joint  kinematic  commands  that  give  rise  to  zero 
Coriolis.  In  this  case  it  is  desired  to  solve  the  following  optimization  problem  by  minimizing 
JBi  with  respect  to  a  scaler  control  variable  p(t)  where: 


r 

JB,  =  (l/2)/[|  I  0(z)  -©k(z)  I  Iri  +  (1  -  p(z))2  r2]  dz 

J  to  (25) 

Here,  again,  the  variable  0,^  becomes  the  state  vector  and  the  control  variable  (or  variable  to 
be  determined)  is  p(t).  The  time  interval  t  in  [  to,  tj  describes  the  time  duration  in  which  the 
solution  is  desired.  If  the  motion  simulator  is  viewed  as  a  robotic  system,  then  ©(.)  is  a 
vector  of  joint  angles  that  are  desired  and  specified  via  equation  (19)  since  W  is  known  a 
priori.  ©^  represents  the  actual  joint  angle  time  histories  resulting  from  joint  motions  with 
reduced  Coriolis.  p(t)  will  be  defined  as  a  result  of  the  linear  quadratic  optimization 
procedure.  Note  R1  is  a  positive  definite  symmetric  matrix,  r2>0,  and  the  norm  denotes 
Euclidian  norm.  Minimizing  JBj  is  subject  to  the  dynamic  constraint  (the  state  vector 
equation  for  ©^  which  is  dependent  on  the  control  variable  p(t)): 

(d/dt)  ©k  =  p(t)  [J2]-'  W  (26) 

where  p(t)  (not  equal  to  1)  is  to  be  selected.  Equation  (26)  can  be  numerically  integrated 
forward  to  determine  the  commanded  or  kinematic  commands  ©k(t).  If  p(t)  =  l,  then  the 
solution  of  (26)  is  given  by  equation  (19)  and  contains  Coriolis.  Note  that  ©ic(t)  is  described 
via  equation  (26)  since  W(t)  is  specified  a  priori.  The  optimization  of  (25)  subject  to  (26) 
also  has  a  state  variable  path  constraint  during  t  in  [t^,  tj  of  the  form: 

(d/dt)W,(t)  =  (d/dt)W,(t)[©k(t),  (d/dt)©k(t)]  =  0  (27) 

It  is  noted  that  the  hard  constraint  on  the  attitude  variable  in  equation  (27)  is  what  gives 
method  1  the  classification  as  zero  Coriolis  in  attitude  space. 


13 


A  second  solution  to  this  problem  is  as  follows: 

Method  2  -  Minimum  Coriolis  in  an  Attitude  Space 

A  second  method  to  approach  this  problem  is  to  determine  a  modified  time  history  ©^(t)  to 
produce  a  command  signal  for  the  three  actuators  of  the  motion  simulator  which  involves 
some  nonzero  Coriolis,  but  reduced  from  the  expression  specified  in  equation  (19).  This 
approach  to  the  problem  differs  from  Method  1  because  the  main  concern  is  the  reduction  of 
Coriolis  precisely  at  the  body  axis  system  as  it  affects  the  human  subject  with  respect  to  the 
orientation  variables.  This  solution  also  allows  for  the  existence  of  some  nonzero  Coriolis, 
which  is  the  case  in  both  the  simulator,  as  well  as  the  actual  aircraft.  In  this  formulation,  the 
objective  is  to  minimize  JB3  with  respect  to  a  variable  p(t)  where: 


JBj  =  (l/2)/[l|(d/dt)W,  [eK(z),(d/dt)ek(z)]||R,  +  (l-p(z))2r2]  dz 
to 


subject  to  the  dynamic  relationship: 
(d/dt)©,  =  p(t)  [J2]-i  W(t) 


(28) 


(29) 


Again,  as  in  Method  1,  the  state  vector  becomes  ©^(t)  and  the  control  variable  is  p(t).  Note 
that  (d/dt)W<,  is  defined  in  equation  (21)  and  thus  known  explicitly.  Hence,  in  Method  1 
there  exists  a  hard  constraint  on  the  Coriolis  terms.  In  Method  2,  however,  the  optimization 
procedure  is  used  to  minimize  these  terms,but  some  artifact  is  permitted  (resulting  in  the 
existence  of  nonzero  Coriolis)  through  the  tradeoff  between  the  attitude  variables  and  the 
joint  commands.  The  appendix  describes  in  detail  the  actual  optimization  procedures  for 
each  of  the  two  solutions  to  these  problems.  These  results  are  summarized  in  the  next 
section. 


Results  of  the  Optimization  Procedure 

With  reference  to  the  appendix,  using  the  assumption  that  the  orientation  Jacobian  J2(©)  was 
nonsingular  and  could  be  used  as  a  simplification  for  the  determination  of  the  p(t)  equation 
(if  not,  J2  (©k)  would  be  a  more  accurate  approximation, or  some  damped  least  squares 
method  could  be  employed),  it  was  also  noted  that  in  both  solutions,  the  expression  for  p(t) 
was  found  to  be  of  the  form: 

p(t)  =  1  +  (Coriolis  reduction  term)  (30) 
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because  if  p=l,  then  full  Coriolis  exists.  The  second  term  in  equation  (30)  must  reduce  the 
Coriolis.  This  second  term  was  found  to  be  (in  the  solution  of  both  optimization  problems): 

Coriolis  reduction  term  =  -(r2)  *  (J2)  ‘  W  (31) 

where  L  is  an  adjoint  variable  which  will  be  eliminated  in  the  ensuing  analysis. 

Thus,  by  judiciously  selecting  r2,  one  can  adjust  the  Coriolis  effects  to  be  made  arbitrarily 
small.  Of  course,  a  tradeoff  would  then  exist  between  minimization  of  the  Coriolis  terms 
compared  to  deviations  of  from  the  desired  0  values. 

To  illustrate  the  applicability  of  this  approach,  the  second  method  will  be  investigated  via  a 
computer  study  of  the  mathematical  model  describing  the  motion  simulator  (centrifuge). 


Simulation  Results 

Applicability  of  The  Approach  to  a  Motion  Simulator: 

To  apply  this  approach  to  a  motion  simulator,  it  is  first  required  to  have  available  the 
orientation  Jacobian  for  the  centrifuge  simulator  with  respect  to  the  body  axis  system  in 
Figure  3.  In  (Repperger,  [1,6]),  The  J2  (Orientation  Jacobian)  was  determined: 

SsyC„,-s,t  cl+c,yC„sl,  slc2s,t+clc2c,yC„+s2s^C3, 

J2=  -c,y,  si  s^,  cl  c2  s,y  -  c.y  s2 

“SjySjt,  ~clCj(“  slCjySgf,  slc2  Cg(~clc2CgySgt  “s2  Sjy  Sj{  (32) 

where  the  subscripts  sy  and  st  again  refer  to  the  seat  yaw  rotation  and  seat  back  angle, 
respectively  as  depicted  in  Figures  1,3.  For  the  example  chosen  in  this  report,  the  choice  is 
made  of  0gy=  Pi/2,  ©5t=©.  Geometrically  this  corresponds  to  the  seat  facing  tangentially 
in  the  direction  of  rotation  of  the  centrifuge. 

For  the  seat  configuration  specified  above,  the  Jacobian  now  becomes: 

1,  0,  s2 

J2  =  0,  si,  cl  c2  (33) 

0,  -cl,  si  c2 

To  illustrate  the  applicability  of  this  approach,  it  will  be  demonstrated  on  a  special  motion 
trajectory  termed  the  "Cobra." 

The  Cobra  is  a  "supermaneuver"  (aircraft  agility  maneuver  usually  involving  high  angle  of 
attack  and  post-stall  motions)  which  has  been  demonstrated  at  airshows  (Aviation  Week, 
[10],  and  Cook,  [11])  and  has  special  significance  in  a  combat  scenario  (cf.  Figure  6). 
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When  the  lower  aircraft  performs  this  supermaneuver  in  Figure  6,  it  pitches  up  rapidly  and 
slows  down  in  the  sky  as  a  tactical  ploy  to  confuse  an  opposing  aircraft  and  to  point  first. 

As  the  lower  aircraft  translates  along  its  velocity  vector  with  its  nose  pitched  up,  it  resembles 
a  cobra  snake  ready  to  strike.  It  is  simulated  with  the  assumptions  of  1000  feet  altitude  and 
a  velocity  change  of  from  300  Miles  Per  Hour  to  150  Miles  Per  Hour.  Most  applications  of 
this  supermaneuver  have  been  performed  at  low  altitude  and  slow  speed  due  to  the  potential 
structural  damage  to  the  aircraft  during  this  flying  procedure.  The  following  commands  on 


the  joint  motions  of  the  centrifuge  will  emulate  this  motion: 

Roll  Rotation  =  ©3(t)  =  (Pi/ 180) (30-2. 5 1)  for  t  in  [0,5] 

=  (Pi/180)*(5+2.5t)  for  t  in  [5,10]  (34) 

Yaw  Rotation  =  0.0 

Pitch  Rotation  =  -02(0  =  (Pi/3)*[l-cos(Pin/5)]  (36) 

Arm  Velocity  =  (d/dt)©i(t)  =  8.8-0.66t  in  RPM  for  t  in  [0,5]  (37a) 

=  5.5  +  0.66(t-5)  in  RPM  for  t  in  [5,10]  (37b) 


Figure  6.  Pougochev’s  "Cobra"  Supermaneuver 

Note  that  the  ©2  and  ©3  variables  have  units  of  radians,  but  (d/dt)  ©1  is  in  RPM.  From 
equations  (37a-b),  it  is  obvious  that  (d/dt)  ©,  being  nonzero,  is  changing  with  time.  Hence, 
this  simulator  has  Coriolis  because  a  frame  (the  fork  =  ©2,  see  Figure  1)  rotates  on  another 
moving  frame  (the  arm  =  ©,)  and  the  cab  =  ©3  rotates  on  a  moving  frame  (the  fork  = 
©2).  To  calculate  the  magnitude  of  the  Coriolis  accelerations,  it  is  necessary  to  mainly 
consider  the  inertial  frame  vector  Jo  term.  This  term  is  the  most  disturbing  to  the  pilot  in  this 
supermaneuver  because  it  is  in  a  direction  (body  axis  Gy)  to  produce  disorientation  effects. 
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Table  3  lists  the  respective  components  of  BB3  as  specified  in  equation  (23)  for  t  in  [0,  10]; 

Table  3.  Relative  Components  of  BB3(t)  for  the  Cobra  Maneuver 

(d/dt)ei=Wl,  (d/dOGj  (d/dt)©3 

tin  [0,5]  2.5-0.24t  -(PiVl5)[sin(Pi  t/5)]  -0.105 

tin  [5,10]  1.3+0.24t  -(PP/15)[sin(Pi  t/5)]  0.105 

Using  equation  (23)  for  BB3,  Figure  7  illustrates  the  Coriolis  acceleration  for  this 
supermaneuver  versus  time  for  t  in  [0,10].  The  units  for  the  dependent  variable  are  Gy  (1 
Gy  =  32.2  ft./s/s).  From  Figure  7  it  is  seen  that  the  pilot  experiences  a  Coriolis  Gy  of  over 
+.  1.0  Gy  during  this  supermaneuver.  This  is  a  sufficiently  large  value  of  motion  artifact,  as 
perceived  by  the  pilot,  as  to  make  the  centrifuge  simulation  unrealistic. 

To  correct  this  situation,  the  joint  commands  will  be  integrated  forward  using  equation  (29) 
and  an  optimization  procedure  is  applied  similar  to  Method  2  with  the  objective  of 
minimizing  Coriolis  in  attitude  space.  For  the  objective  function  JBj  in  equation  (28),  the 
choice  is  made  of  R1  =  21,  r2  =  8,  1,  and  0.2.  This  results  in  the  following  two  point 
boundary  value  problem: 


L(tf)  =  0  (38) 

(d/dt)L’'  =  -(1/2)  *  (d/d  ek)[(d/dt)W7  R1  (d/dt)WJ  (39) 

p(t)  =  1  -  (r2)-‘  (J2)->  W  (40) 

(d/dt)©k  =  p  (J2)-‘  W(t)  (41) 


BB3  Coriolis  For  The  Cobra  Maneuver 

Centrifuge  Induced  Coriolis 


- Centrifuge  Induced 

Figure  7.  B  Bar  3  Coriolis  for  the  Cobra  Supermaneuver 


17 


It  is  noted  that  the  term  -(l/2)*(d/d  ©k)[(d/dt)W7  R1  (d/dt)WJ  is  obtained  from  equation 
(23)  for  BB3  =  (d/dt)  We.  This  assumption  is  made  because  the  main  concern  is  to 
minimize  accelerations  in  the  lateral  Gy  direction  for  each  of  the  L(t)  components.  The  Gy 
component  is  the  most  disturbing  to  humans  and  it  is  desired  here  to  change  the  influencing 
functions  L(t)  as  per  this  Gy  component.  If  the  interest  was  to  reduce  Coriolis  effects  in  the 
other  two  axes,  it  would  be  a  simple  process  to  include  the  terms  AB3  and  CB3  in  the 
ensuing  analysis.  The  equations  for  the  L(t)  variables  now  become  (with  the  assumption  that 
R1  is  diagonal  with  diagonal  elements  rl): 


(d/dt)L,  =  -(1/2)  *  ((d/d  ©Ik)  rl  [BB3f) 

(42) 

(d/dt)L2  =  -(1/2)  *  (d/d  ©aJ  rl  [BB3]' 

(43) 

(d/dt)  Lj  =  0 

(44) 

m  =  0 

(45) 

The  differential  equation  for  Lj  has  zero  on  the  right  hand  side  due  to  the  independence  of 
AB3,  BB3,  and  CB3  on  ©3.  Since  the  vector  L(tf)  =  0,  this  leaves  L,  and  Lj  as  the  only 
nonzero  variables.  The  relationship  for  BB3  is  substituted  from  equation  (23)  and  the  partial 
differentiation  is  first  taken  and  then  the  adjoint  equations  are  then  integrated  backwards. 
Since: 


(d/dt)  ©k  =  p(t)  [J2]-‘  W(t)  (46) 

and: 

p(t)  =  1-  (r2)-‘  (J2)-'  W  (47) 

This  corrects  the  ©;  terms.  The  ©i  are  then  integrated  one  step  forward  to  obtain  their 
respective  corrections.  The  Coriolis  can  then  be  recalculated  with  these  changes  in  the  joint 
commands. 

Figure  8  illustrates  1^(0,  the  adjoint  variable  plotted  for  the  case  Rl  =21.  Figure  9  is  a  plot 
of  the  optimal  control  gain  p(t)  for  the  case  r2=0.2.  In  Figure  10,  plots  of  the  term  BB3  are 
illustrated  for  the  choice  of  the  original  Coriolis  and  r2=l  which  was  the  optimum  weighting 
selected  from  this  numerical  study.  From  Figure  10  it  is  seen  that  a  substantial  reduction  of 
Coriolis  is  possible.  Comparing  over  time  the  sum  of  the  squares  of  the  BB3  terms,  the 
minimum  Coriolis  is  calculated  for  r2=1.0  and  compared  to  the  original  Coriolis  (resulting 
in  a  65.4  percent  reduction  in  Coriolis  accelerations  in  the  Gy  direction  produced  by  the 
motion  simulator  when  it  adapts  these  inverse  kinematics  algorithms). 
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If  the  Coriolis  component  is  reduced,  there  exists  a  tradeoff  in  the  adapted  kinematic 
commands  from  the  desired  commanded  joint  time  history  ©(t).  Examining  equation 
(28)  with  the  two  weighting  matrices  R1  and  r2  and  utilizing  the  fact  that  the  vector  ©(t)  is 
3x1,  the  cost  (in  a  norm  sense)  of  reducing  this  Coriolis  component  was  determined  to 
produce  a  compromise  in  the  joint  fidelity  commands  of  only  5.1  percent.  This  amount  of 
loss  of  motion  fidelity  may  not  be  perceived  by  humans  in  this  complex  acceleration 
environment. 


The  Adjoint  Variable  Versus  t 


—  Adjoint  Variable 


Figure  8.  The  Adjoint  Variable  Versus  Time  -  Lambda  2  (t) 


p{t)  =  The  Optimal  Gain 

Versus  Time 


- p(t)  Versus  Time 

Figure  9.  p(t)  -  The  Optimal  Gain  Versus  Time 
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BBS  Coriolis  For  The  Cobra  Maneuver 

The  Effect  of  r2  on  the  Solution 


- B3  Original  - B3  (r2=1) 

Figure  10.  B  Bar  3  Coriolis  for  the  Cobra  Supermaneuver  -  The  Effects  of  r2  on  the 
Solution 


SUMMARY  AND  CONCLUSIONS 

A  three  axis  motion  simulator  is  modeled  as  a  robotic  manipulator.  The  equations  of  motion 
are  derived  and  the  terms  due  to  Coriolis  are  quantified.  Two  optimization  algorithms  are 
presented  to  minimize  the  Coriolis  produced  in  the  simulator  at  the  end  effector.  It  is  shown 
with  minimal  reduction  of  motion  fidelity,  that  over  65  %  of  Coriolis  can  be  reduced  with 
proper  input  commands  to  modify  the  motion  field  produced.  When  this  reduction  of 
Coriolis  is  realized,  the  motion  fidelity  is  only  compromised  about  5%  which  may  not  be 
perceivable  by  the  human  subjects  that  participate  in  the  human  runs  in  the  mechanical 
system  considered  herein. 
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APPENDIX 


Method  1  falls  within  the  context  of  a  linear  quadratic  optimization  problem  with  path 
constraints  along  the  optimal  trajectory  (Bryson  and  Ho,  [9]).  Method  2  is  also  a  quadratic 
optimization  problems,  with  the  control  entering  in  a  linear  manner.  We  now  reformulate 
each  optimization  problem,  defining  the  state  variables,  the  control,  and  the  procedure  to 
solve  each  problem. 

Method  1  -  Minimize  with  respect  to  the  control  variable  p(t): 


JB,  =  (l/2)ri^|  ©(z)  -e,(z)  I  +  (1  -  p(z)f  r2  dz  (48) 

Jto 

where  the  state  variable  dynamic  constraint  is  specified  by; 

(d/dt)  ©k  =  p(t)  [J2]-‘  W  (49) 

and  state  variable  path  constraint  is  given  by: 

C  =  (d/dt)W,(t)  =  (d/dt)W,(t)[©k(t),(d/dt)©,(t)]  =  0  (50) 

The  solution  of  this  problem  is  specified  once  the  Hamilitonian  H(.)  is  chosen.  For  this 
case: 


H  =  1/2  *  ((©-©k)^  R1  (©-Ok)’")  + 

(1/2)*  (1-p)^  r2  +  [  p  (J2-')]  W  +  nu  C  (51) 

where  the  L(t)  vector  is  the  adjoint  state  variable  and  nu(t)  is  another  influencing  function 
chosen  to  deal  with  the  equality  constraint  in  equation  (50).  It  is  noted  that  in  equation 
(49)  the  control  p  is  linear  in  the  state  equation,  and  that  W(t)  and  0(t)  are  known  and 
specified  a  priori.  The  necessary  conditions  for  optimality  require; 

(d/dt)  L''(t)  =  -  Hek  (52) 

with  transversality  condition; 

L(tf)  =  0  (53) 

Thus: 

(d/dt)  V  =  -  (©-©k)^  R1  -  nu  (d/d  ©k)  [C]  (54) 
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The  second  necessary  condition  is: 

Hp  =  0 


(55) 


which  yields: 

(p-1)  r2  +  L’’  (J2)-'  W  =  0  (56) 

which  results  in: 

p(t)  =  1  -  (r2)-‘  (J2)-‘  W  (57) 

Thus,  a  two  point  boundary  value  problem  is  established  in  which  the  adjoint  equations  for 
L(t)  are  integrated  backwards  using  equations  (53)  and  (54)  with  p(t)  is  calculated  forward 
from  equation  (57),  Interestingly  enough:  p=l+  (a  second  term).  If  p=l,  then  Coriolis  is 
present.  The  second  term  must  be  the  Coriolis  reduction  term  which  adds  linearly.  The 
solution  for  problem  2  is  now  presented  because  it  minimizes  the  Coriolis  in  the  motion  field 
directly  experienced  by  the  human  subject.  It  also  allows  for  some  nonzero  Coriolis  artifact 
which  is  the  case  for  the  simulator,  as  well  as  in  the  actual  aircraft. 

Method  2  -  Minimize  with  respect  to  the  control  variable  p(t): 

r 

JB3  =  1/2/1  |d/dt(wj[e,(s),(d/ds)  e,(s)]|  Ir,  +  (1-  p(s)y  ri  ds 

/  (58) 

subject  to  a  dynamic  constraint  on  the  state  variables  of  the  form: 

(d/dt)  ©k  =  p(t)  [J2]-‘  W  (59) 

For  this  case  a  Hamilitonian  is  defined  as  follows: 

H  =  1/2  (d/dt)W/  R1  (d/dt  (We))  +  1/2  (l-p)^r2  +  L‘^[p(J2)-‘W]  (60) 

The  necessary  conditions  for  optimality  require: 

(d/dt)L’-(t)  =  -  H(ek)  (61) 
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with  transversality  condition: 

L(h)  =  0 


(62) 


Thus: 


(d/dt)  =  -  1/2  (d/d©J[(d/dt)  (W,)^  R1  (d/dt)  WJ  (63) 

The  second  necessary  condition  (Hp  =  0)  yields; 

(p-1)  r2  +  (J2)-*  W  =  0  (64) 

hence  the  familiar  control  law  results  (identical  to  equation  (57)): 

p(t)  =  1  -  (r2)-‘  (J2)-‘  W  (65) 

since  (d/dt) Wc(©k)  is  known  (cf,  equations  (21-24)),  this  problem  is  easily  solved. 

Finally  it  is  noted  that  it  has  been  assumed  that  (J2[©])  '  is  approximately  (J2[©J)  ‘  to 
eliminate  the  complexity  of  this  factor  into  the  analysis.  If  the  (d/d  ©k)  (J2)'*  term  were  to 
appear  in  any  of  the  L(t)  equations,  the  problem  becomes  increasingly  more  difficult. 
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